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ABSTRACT 

We consider a simple gravitational-heating mechanism for the long-term quenching 
of cooling flows and star formation in massive dark-matter haloes hosting elliptical 
galaxies and clusters. The virial shock heating in hal oes ^ 10^^ A/0 triggers natural 
quenching in 10^^~^^Mq haloes ([Birnboim et al.ir2007t ). Analytic estimates and simple 
simulations argue that the long-term quenching in haloes ^ Mmin '--^ 7 x lO^^M© could 
be due to the gravitational energy of cosmological accretion delivered to the inner- 
halo hot gas by cold flows via ram-pressure drag and local shocks. Mmin is obtained 
by comparing the gravitational power of infall into the potential well with the overall 
radiative cooling rate. The heating wins if the gas inner density cusp is not steeper 
than r~^^^ and if the masses in the cold and hot phases are comparable. The effect 
is stronger at higher redshifts, making the maintenance easier also at later times. 
Particular energy carriers into the halo core are cold gas clumps of ^ 10^~*Mq. 
Clumps ^ IO^Mq penetrate to the inner halo with sufhcicnt kinetic energy before 
they disintegrate, but they have to be ^ 10^ Mq for the drag to do enough work in a 
Hubble time. Pressure confined 10*K clumps are stable against their own gravity 
and remain gaseous once below the Bonnor-Ebert mass ~ IO^Mq. Such clumps are 
also immune to tidal disruption. Clumps in the desired mass range could emerge by 
thermal instability in the outer halo or in the filaments that feed it if the conductivity is 
not too high. Alternatively, such clumps may be embedded in dark-matter subhaloes if 
the ionizing flux is ineffective, but they separate from their subhaloes by ram pressure 
before entering the inner halo. Heating by dynamical friction becomes dominant for 
massive satellites, which can contribute up to one third of the total gravitational 
heating. We conclude that gravitational heating by cosmological accretion is a viable 
alternative to AGN feedback as a long-term quenching mechanism. 
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1 INTRODUCTION 

The observed properties of red & dead elliptical galaxies 
indicate a robust quenching of star formation above a 
threshold corresponding to halo mass ~ IO^^Mq, effec- 
tive especi ally after a characteristic redshift z ~ 2 (ref- 
erences in lDekel fc Birnboimll2006l ). Similar quenching 
is required in order to suppress cooling flows in clus- 
ters, i.e., in halos ~ IO^-^^Mq fPabian 1994). In turn, 
semi-analytic simulations show that the introduction of 
quenching above a threshold halo mass (sometimes ac- 
companied by a critical black-hole or bulge mass) is 
a crucial element in any model that tries to fit the 
data concerning the distribution of galaxy properties 



(ICroton et al.ll2006t ICattaneo et alll2006t iBower et al.l 
120061) . 

The desired quenching mechanism should provide 
the required energy over time and explain its coupling 
to the whole gas reservoir in the inner dark halo. It 
needs to address both the trigger of quenching and its 
long-term maintenance for many Hubble times, as well 
as the characteristic mass and redshift associated with 
it. 

There is an ongoing intensive effort to study "AGN 
feedback" as the potential source of quenching. At 
this point, the physics of AGN feedback, e.g., how 
it couples to the extended halo gas and how it pro- 
vides long-term maintenance, are still difficult open 
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issues (e.g. iBriiggen et al.l [20071 ). The bright quasars 
have been argued to be capable of providing long- 
term feedback via radiation despite their short life- 
times (jCiotti fc Ostrikeill2007D . but perhaps more nat- 
ural sources for self-regulated AGN feedback are the 
weaker radio-loud A GNs that r adiate for a long time at 
low power (review by'Bcst 2007, and references therein). 
In either case, the desired characteristic mass and red- 
shift do not emerge in an obvious way from the black- 
hole physics. 

Gaseous major mergers, which were suggested as 
the trig ger for quenching vi a starbursts or quasar acti- 
vation ( Hopkins et al.l [20071 ). also have hard time pro- 
viding a satisfactory explanation for the characteristic 
mass, and it is not clear that their frequency and the as- 
sociated s tarburst efficienc ies are sufficient for the pur- 
pose I'e.g. lLotz et all[200l ?: lCoxet al.l[200l . 

Here we pursue a preliminary feasibility test of a 
very simple alternative, that the gravitational energy 
associated with the cosmological accretion of baryons 
into the dark-matter haloes is the dominant source of 
quenching, responsible both for its trigger and its main- 
tenance. The trivial fact that the baryons tend to segre- 
gate from the dark matter and condense at the bottom 
of halo potential wells implies that they bring in energy 
that must dissipate and coul d keep the halo gas hot. 

Our earlier analy sis (jBirnboim fc Dekel l2003t 
iDekel &: Birnboim[ 120061 ). based on either analytic cal- 
culations or simulations in the idealized case of spherical 
symmetry, explained the existence of a robust thresh- 
old halo mass for a stable shock at the virial radius, 
Mihock ^ 1O^^M0. This is a refinement of the classical 
idea attributing galaxy formation to cooli ng on a dy- 
narnical t ime scale (Rees & Os trikcr 1977: B i nneylll977l : 
Siigi 19771: [Wh ite & Rccs 19781 : iBlumenthal et al.lll984t ). 



Birnboim fc Dekel (,2003i) showed that in halos below 



the threshold mass, rapid cooling prevents the post- 
shock gas pressure from supporting the shock against 
gravitational collapse. In this case, the gas flows cold 
(~ lO'^K) into the inner halo, where it may eventu- 
ally shock, form a disk and possibly generate an ef- 
ficient mode of star formation. When the halo grows 
above Mghock, a stable shock propagates toward the 
virial radius, halting the inf ailing gas and creating a hot 
medium in quasi-static equilibrium at the halo virial 
temperature. The transition to shock stability occurs 
when the standard radiative cooling time equals the 
time for compression behind the shock, tcompross = 
{21/5)p/p ~ (4/3)i?/F, where p is the gas density be- 
hind the shock, R is the shock radius and V is the infall 
velocity into the shock. Early indications for cold accre- 
tion versus shock virialization as a function of redshift 
origin ated from three-dimensional cosmological simula - 
tions (jKatj 119921: iKav et al.|[2000l : iFardal et al.ll200lD . 
Advanced cosmological simulations solidified the tran- 
sition to a hot medium at a critical halo mass in par- 
allel with the theoretical understanding of the effect 
(iKatz et al.l [20031: iKeres et al.ll2005t iDekel fc Birnboim! 
[20061 : [Birnboim et al. 2007 . fig. 1, based on simulations 
by A. Kravtsov). This has become the standard lore 
used in modeling of gala.xy formation ([Croton et al.[ 
[20061 [Cattaneo et al.[[2006[ : [Bower et"al[[2006t) . 



The introduction of a hot medium filling the halo 
is a natural trigger for quenching of cold gas supply. 
This is either directly, by shock heating of all the gas, 
or indirectly, by providing a hot dilute medium that is 
vulnerable to heating by sources such as AGN feedback. 
This very naturally explains the origin of the charac- 
teristic thre shold halo mass re sponsible for the galaxy 
bimodality pekel fc Birnboim|[2006 ). 

While the role of virial shock heating as the trigger 
for quenching and as a necessary condition for mainte- 
nance arc becoming widely appreciated, we address here 
the novel idea that the accretion itself may also be the 
dominant source of energy for long-term maintenance, 
without appealing to AGN feedback or in parallel to it. 

Using spherical hydro simulations, Birn boim et al.[ 

||2007( ) have found that a uniform virial accretion at 
the average cosmological rate leads by itself to efi'ec- 
tive long-term quenching starting at z ^ 1 in halos 
of masses 2 x 10^^ M < IO^^Mq today, namely 
in groups of galaxies. This is a generic sequence of 
events, termed SAMBA (for Shocked- Accretion Massive 
Burst and Shutdown) , due to a moderate accretion rate 
through a rapidly expanding virial shock. The long-term 
quenching follows an earlier tentative quenching phase, 
when the halo is ^ lO^^M©, and a subsequent, rapid, 
massive accretion of ~ 10^^ Mq of gas into th e inner 
halo . Cosmological siniulation s by [Naab et all (|2007f) 
and iLibeskind fc Dekel (|2007[ ) provide first clues that 
such a mechanism may actually work once the numer- 
ical r esolution is suffic i ent for resolving s mall clumps 
(also [Motl et all [20041 : [Burns et all [2007t [Nagai et al.[ 
[2001 . 

However, [Birnboim et al.l (|2007f ) found that uni- 
form spherical accretion does not provide long-term 
quenching in halos more massive than 10^'^ Mq. This 
is because the accretion rate was too high at the early 
time when the halo was ~ lO^^Af©. We address here 
a more robust phenomenon that is also valid in more 
massive haloes and in realistic cosmological geome- 
tries, where the gravitational energy of accretion is the 
poten tial source of long-ter m qu enching maintenance. 
Both [Birnboim et all (|2007l ) and [Khochfar fc Ostriked 
([20071 ) have made preliminarv attempts at addressing 
the potential of a scenario along these lines, assuming 
that the energy is transferred to the inner-halo gas via 
clumps. It is important to note that heating by dy- 
namical friction acting on subhalocs, which is a non- 
negligible energy source, is in sufficient for balancing 
the cooling rate ([El-Zant et all ,2004: Kim et al.l[2005[ : 



iSiiacki fc Soringell 120061: iKhochfar fc Ostriked \20m . 
We rather focus here on the effect of ram-pressure drag 
(or shocks) exerted by a hot ambient gas on cold gas 
clumps 13 

In ^ we find analytically that the overall gravita- 
tional accretion power is adequate for keeping the halo 



^ We refer in this paper to ~ lO^K clumps as "cold", to dis- 
tinguish them from "warm" gas at ~ 10^ K and "hot", X-ray 
emitting gas of ^ lO^K. Similar clumps are sometimes referred 
to as "warm", reserving the term "cold" to gas that has cooled 
by molecular cooling to lO^K. 
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gas hot in halos above a threshold mass ~ 6 x lO^^M©, 
provided that the energy is deposited in the inner halo. 
We evaluate the scaling with redshift, density profiles, 
gas fractions and metallicity. In we work out con- 
straints on the clump masses, ttt-c ^ 10^~®Mq, neces- 
sary for effective heating by drag of the hot gas in the 
inner halo. In 21 present simple simulations that 
confirm the analytic estimates and demonstrate the po- 
tential of keeping the halo gas hot by clumpy accretion. 
After sketching the method of simulation, we analyze 
the time evolution of individual clumps, and use the 
results for the heating-to-cooling ratio inside the halo 
to obtain the constraints on halo and clump masses. 
In ^jni we discuss the possible role of dynamical friction, 
the physics of dissipation via turbulence, and the pos- 
sible origin of clumps by thermal instability of inside 
subhaloes. In ^we conclude our findings. 



2 GLOBAL HEATING VERSUS COOLING: 
MINIMUM HALO MASS 

The gravitational energy associated with infall into 
massive haloes is clearly a viable source for quench- 
ing because, by virtue of the virial theorem itself, it 
is comparable to t he overall ther mal energy (see also 
IWang fc Ab"eill2007[ ). As argued bv lKhochfar fc Ostriked 
( 2007 ). this gravitational energy can compete with 
the popular energy source of feedback from a central 
AGN. For a total stellar mass M*, the gravitational 
energy is ^ 0.5M^V^^^, while the AGN feedback can 
be expressed as ^ eM^c^, so the ratio of the two is 
- (V;sc/400kms-i)2(e/10-6)-i. This is of order unity 
for typical escape velocities in big ellipticals and clus- 
ters and for typical AGN efficiencies. Unlike the AGN 
energy which is likely to heat the densest gas and largely 
get radiated away, the gravitational heating can be de- 
posited much more uniformly so less energy is wasted 
as radiation. 

A necessary condition for the accretion to pro- 
vide an effective quenching mechanism is that there is 
enough gravitational power in the accretion to overcome 
the overall rate of energy loss by radiative cooling. By 
considering first this global requirement we will obtain 
here a robust lower limit for the halo mass in which 
gravitational quenching is possible. In the following sec- 
tions we will obtain limits on the clump masses by con- 
sidering the energy balance in the inner halo, with the 
heating provided specifically by ram-pressure drag on 
cold clumps. 



2.1 Halo Equilibrium Profiles 

In order to evaluate the rates of heating versus cool- 
ing we use a simple model for gas in hydrostatic equi- 
librium within the potential well of a spherical dark- 
matter halo. The halo virial mass Mv and radius i?v 
are defined in the standard way such that the mean 
mass over-density within i?v relative to the cosmologi- 
cal background is A '-^ 200. Recall that the cosmological 
density is p ~ 2.76 x 10~^°Qaio.3hl a~^gcm~^ , with 
a = {1 + z)~^ the universal expansion factor, and in a 



flat univers e A varies from — 180 a t 2: > 1 to :^ 340 at 
z = (e.g. iBrvan &: Normanlll998t iDekel fc BirnboimI 



I2OO6I Appendix A). The mean mass density within the 
virial radius is then expressed conveniently as 

Pv 



5.52 X lO^^M^^gcm"^ , 



(1) 

(2) 



(4^/3)i?3 

A = (A200 ho.ry^^^ ' 

With the virial velocity = GA'U/Rv, writing M13 = 
Mv/IQI^Mq, Ri = i?v/Mpc and V300 = K/BOOkms^i, 
the standard virial relations can be written as 

Mi3 ~ 3.42 RlA-^ ~ 1.64 Vio^A^/'^ . (3) 

We assume that the total and the hot-gas density 
profiles are both of a generalized NFW functional form, 

Ps _ r 



p{r 



(4) 



x^il + x)^-"' rs' 

They can be parameterized by the inner slope a, the 
mass M within the virial radius, and the concentration 
parameter C = Rv/r^. We denote the total hot-gas mass 
within the virial radius by Mg = fgMv The concentra- 
tion, for each of the two components, is assumed to 
follow the cosmolog ical average for dark-matter haloes 
(iBullock et al.ll200lb . 

(5) 



C~ 9.0M 



-0.15 



13 



Hydrostatic equilibrium requires at every radius r 



GM{r) 



1 



d 



Pg{r) dr 



k 



Pg{r)T{r) 



(6) 



where M{r) is the total mass within r, Pg{r) and T(r) 
are the gas density and temperature profiles, k is the 
Boltzmann constant and m ~ 0.59TOp with mp the pro- 
ton mass. Thus, for given density profiles, we can eval- 
uate the temperature profile via 

M{r) 
^rpg{r) — ^ 



T(r) 



Gm 1 



k pg{r) 



TV, (7) 



with the external boundary condition set by the virial 
temperature at i?v 



m 2 
or in the cosmological context 



Te ~ 3.23V4 



2.33 M^l^A~ 



T 

lO^K 



The entropy profile is then proportional to 
T{r) 



K{r) cx 



Pg(r)2/3 



(8) 



(9) 



(10) 



In our fiducial model we assume an NFW cusp for 
the total mass, ad = 1, and a constant-density core for 
the gas, ckg = 0, with /g = 0.05. The resultant pro- 
files are shown in Fig. [1] The gas temperature is only 
weakly varying about the virial temperature. The en- 
tropy is rising at large radii and is approaching a con- 
stant entropy floor at small radii, in qualitative agree- 
ment with t he typical profiles indicated by X-ray ob- 
servations (jDonahue et al.l [20061 : iPratt et al.ll2006() and 
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Figure 1. Halo profiles in hydrostatic equilibrium. Total density 
(blue), gas density (magenta), gas temperature (long-dash red) 
and gas entropy (short-dash green) as a function of radius, all 
in terms of the virial quantities. The total mass is IQl^M© at 
2 = 0. The total density profile is NFW. The hot gas fraction is 
/g = 0.05. The fiducial case is with a gas core and entropy floor 
(top), compared to a case with an NFW gas cusp (bottom). 



with tho se obtained in hydrodyn amical simulations of 
clusters ([Faltenbacher et ani2007f ). Shown for compar- 
ison is a case with a gas cusp, ag = 1, for which the 
temperature is still roughly constant, but the entropy 
is improperly rising even at small radii. 

For the purpose of analytic estimates, we note that 
the generalized NFW mass profile of eq. ^ can be writ- 
ten as 



M (r) = An ^„ (x) , Ps ^ Pv 



where for an a = 1 cusp 
Ai{x) = \n{x + 1) - 



x + 1 



(11) 



(12) 



Ao{x) = ln(a:; + 1) 



x{3x + 2) 
2(x+ 1)2 



(13) 



2.2 Heating Rate by Gravitational Accretion 

We write the available gravitational accretion power as 



heat 



|A'/'(i?v,r) 



1 



1, 



Vj Me 



(14) 



The first term represents the gain in potential energy 
by infall of clumpy gas from the virial radius to an inner 
radius r, given an accretion rate Mc- The second term 
is the contribution of kinetic energy associated with the 
typical velocities Uj of the clumps at the virial radius. 
The third term represents the energy wasted for un- 
binding infalling cold clumps and heating them to the 
halo virial temperature. Given the density profiles, the 
potential gain is 



|A0| = C^Vl : 



-1 



^^('^'l'dr',(15) 



with V^{r) = GM{r)/r. 

For our fiducial NFW potential well, a^x — l, 



C 



ln(l + x) ln(l + C) 



C 



(16) 



With C = 9, we have Ai ~ 1.49, so (i>{r^) ~ 2.8 and 
0(0) ~ 4.8. In comparison, a singular isothermal sphere 
has 4){r) — ln(i?/r), namely ^(rg) ~ 2.2. 

The average virial accretion rate onto haloes of 
mass Mv at redshift z can be estima t ed us ing the 
EPS formalism, following iNeistein et al. (|200(: ) and a s 
described in Appendix A of 'Birnboi m et al. (|2007f ). 
We use the practical approximation for haloes near 



'1012-i3a/q in ACDM, 



M 



(Mv, t) ~ 0.056 Gyr"' Mi's'^'a 



0.15„-2.25 



(17) 



This is consistent with the estimates from cos- 
mological N-body sim ulations (jWechsler et al.l 120021 : 
INeistein fc Dekel l2007l based on the millennium sim- 
ulation). The accretion rate of clumps is assumed to be 

Thus, for an NFW potential, the heating rate from 
eq. (HI]), in terms of the total thermal energy, is 



heat 



MgK,72 



0.11. 



/c , .r0.15 „-2.25 
7~-'"13 " 
/g 



Gyr- 



(18) 



The strong redshift dependence is due to the higher 
specific accretion rate at higher redshifts, and note the 
relative insensitivity to M^. 

2.3 Radiative Cooling Rate 

The radiative cooling rate per unit mass, at a position 
where the gas density, temperature and metallicity are 
Pg, Tg and Z, is 



and for an a = core 



X 



HTg, Z) pg , 



(19) 
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where x — 0-52 and A(rg, Z) is the atomic coohng func- 
tion. We use the stand ard co ohng function computed by 
[Sutherland fc Dopitai (|1993( ). which is consistent with 
the eq uihbrium state computed by I Gnat fc Sternberg 
(|2007f ). At a radius r within the halo, the cooling rate 
becomes 

(7(r)~7.7xl0-5A_23W/5g(r)/.o5^~'ergs-ig'\ (20) 

where /.05 = /g/0.05 and the gas density profile is ex- 
pressed in a dimensionless way by 



- / N Pg{r) _ J. _ 



(21) 



with pv from eq. ([T]). We define A_23 = 
A/lQ-^^erg cm^s"^ and note that for IO^^Mq, 
the value of A_23 is between 0.6 and 6 for metallicities 
ranging from zero to solar abundance. 

For our fiducial gas profile, ag = 0, we have 



1 



3Ao{C) (l + x) 



3 ■ 



(22) 



With C 9, it is pg(i?v) ^ 0.24, /5g(rs) ~ 30 and 
Pg(0) ~ 243. 

The overall cooling rate of the hot gas is 

£^cooi = 47r/ q{r)pg{r)r^dr. (23) 

Because of the proportionality to p^, the integral is 
dominated by the cooling near the inner radius Vg. 

For our fiducial gas profile, with C ^ 1, and ap- 
proximating Tg{r) = Tv based on Fig. [T] such that 
A is independent of r, we obtain [to an accuracy of 
1-O(10/C3)] 

C3 



qMg 



90Ao(C) 



(24) 



with q the cooling rate of eq. where pg = 1. This 
gives a total cooling rate relative to the thermal energy 



E, 



cool 



0.061- 



r3 



-/.05Mi3'^'A_23^" 



MgT/V2 ^Ao{C)'-"'^'^^'' '^^^ ''^^^^ 

where Cg = C/9 and Ao[C) is given in eq. (fT3|) . e.g., 
Ao(9) ~ 1.00. Note the strong decrease of this specific 
cooling rate with increasing and increasing redshift, 
the latter being dominated by C cx a. 

For comparison with X-ray observations, we note 
that the total cooling rate from our fiducial halo model 
is 



E, 



cool 6.2 X 10^1 fl^ Afi3 A_23^"^ ergs-i 



(26) 



For Mv ^ IQ^'^Mq at z = the virial temperature is 
Tv - 2 X lO^K. Then with Z - 0.3 we have A_23 ~ 3 
so -Ecooi ^ 2 X 10*^ erg s^^. This is indeed compara- 
ble to the observed X-ray output of L^^ ~ 10'*^ergs~^ 
from giant ellipticals and grou ps of virial halo mass 
in the ball-park of lO^-^MfD (iMathews fc Brighenti 
2003"; ' Helsdon fc Ponmanll2003t iHumphrev et al.n200Q : 
Fukaz awaetal J12006). In pa rticular, the observed scal - 
ing with dispersion velocity (jOsmond fc Ponmanll200^ 
can be translated to La, ~ 2 x lO''^ M^g® ergs""'^, quite 



similar to the mass dependence inferred from eq. 
This indicates that our fiducial halo model is roughly 
compatible with X-ray observations, and that in fact 
we have no much freedom in choosing a hot gas frac- 
tion that is very different from the 5% assumed. 

For an analytic estimate of the cooling rate 
one could use the practical approximation by 
iDekel fc Birnboiriil (|2006D to the cooling function. 



A. 



-23 



e.ozo-jTg- 



0.2 



1/2 



(27) 



where Tg is related to the virial mass in eq. ([9|), and 
Z0.3 = {Z/0.3) Zq. We also adopt in our fiducial case 
the redshift evolution of the mean metallicity as crudely 
estimated by Dekel fc Birnboim (2006), 



log(Z/Zo) 



-sz, 



0.3 Zp 



: 



0.17. 



(28) 



2.4 Global H/Q: Minimum Halo Mass 

The ratio of global heating to cooling rate, H/Q, is com- 
puted numerically for any given choice of ofg and a^. 

The heating rate Eheat, as computed from eq. (fT4|) and 
eq. (|15p with Vi — Vv, is divided by -Ecooi, as evaluated 
from eq. ((23)l with eq. ((20|) . Figure [2] shows the resultant 
global H/Q as a function of halo mass. In the fiducial 
case we obtain H/ Q > 1 for 



Mv > M„ 



6.5 X lO^^Mp 



(29) 



This is the lower limit for haloes that permit gravi- 
tational quenching. We see that at Mv > lO^^M©, 
where the natural quenching by shocke d accretion is 
no longer effective ([Birnboim et al.ll2007t ). the quench- 
ing can in principle be provided by gravitational heat- 
ing, once the energy is deposited in the inner halo. At 
Mv = IO^^Mq, H/Q ~ 2.0, while in cluster haloes of 
Mv ~ 1O^^M0, the energy available for gravitational 
heating overwhelms the cooling rate, H/Q ^ 100. The 
almost linear dependence of H/Q on is driven by 
the dependence of the heating rate per unit mass on the 

potential well, c>c Vv^ oc Mv^^, while the cooling rate per 
unit mass is a weaker function of halo mass. 

For the fiducial profiles with ag = and Od = 1, 
we obtain an analytic expression for H/Q by dividing 
eq. (UHl) and eq. ([SS]) . 



H 



CI 



23 ^2.25 



(30) 



We can express H/Q in terms of the parameters Mv, a, 
/c//g, and Zq, by inserting 0(O.li?v) = 3 from eq. (fT6|) . 
C(Mv,a) from eq. (0, /c -^- /g = 0.1, A(Tg,Z) from 
eq. dill with Tg = rv(Mv, A) from eq. Q and Z(Zo, a) 
from eq. (j28p . Using these approximations, we recover 
the fiducial critical value of Mmin — 6.5 x lO^^M©, and 
obtain its crude scaling with the parameters: 

Afniin OC a-^-'', Afmin OC /g//c , ^/min OC Zq"^ . (31) 

Some of these scalings are refiected in the two additional 
curves of Fig. [21 referring to variations of a and fc/ fg 
about the fiducial case. At z = 2, the heating is more ef- 
ficient and we learn that Mmin is reduced to ~ lO^^Mf;^. 
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Figure 2. Estimates of the global gravitational heating rate ver- 
sus cooling rate, H/Q, as a function of halo mass. The fiducial 
case (solid blue) assumes an NFW halo and fg = 0.05 hot gas 
with an inner core {Fig.[T] top). A fraction /c = 0.05 of the mass, 
accreted at the average cosmological rate, is assumed to pene- 
trate to 0.1 -Rv. The metallicity is Z = 0.3^0 and the redshift is 
2 = 0. Shown for comparison are two other cases (dash green): at 
2 = 2, and with fg = 0.075 and /c = 0.025. We see that the over- 
all gravitational heating can in principle overcome the cooling in 
halos ^ 10'^^ Mq. Only haloes abo ve lO'^^M^ (solid red) ha ve the 
necessary shock-heated medium llDekel fc BirnboimI lioOfih ■ and 
haloes in the range lO^'^ — IO^^Mq (dash ed magenta) are natu - 
rally quenched by virial-shocked accretion llBirnboim et al .120071) . 



With a higher fraction of hot gas and a lower fraction of 
mass in clumps, /g — 0.075 and /c — 0.025, we see that 
Mmin — 2 X IO^^Mq, such that gravitational heating 
can barely do the job. 

Figure [3] shows the redshift dependence of Mmin, 
defined by H/Q — 1, for otherwise the fiducial choice 
of parameters. By z — 2, the minimum mass drops 
to ~ lO^^M©, allowing gravitational quenching over 
the whole mass range where a shock-heated medium is 
present. The increase in heating rate due to the higher 
accretion rate (oc a~^'^^) and higher Vv (oc a~^) at 
higher z is almost balanced by the increase in cooling 
rate due to the higher density (oc a~^), leading to the 
apparent weak explicit dependence of H/Q on redshift 
in eq. pO]) . The residual z dependence of H/Q at a 
given Mv is dominated by the decrease of C (cx a) with 
redshift, aided by the corresponding decrease of A via 
Tg and Z. 

In a more realistic dynamical calculation, one can 
expect the over-heating at high z to induce an expan- 
sion of the gas core, reducing /g and Og there. This is 
expected to improve the quenching efficiency at later 
times and thus reduce Mmin there (work in progress) . 

Figure S] displays the dependence of Mmin on the 
inner density profiles of gas and dark matter in equi- 
librium. The dependence on the total density slope 
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Figure 3. Estimated minimum halo mass for gravitational 
quenching, A/,nin, defined by H/Q = 1, as a function of redshift 
2, for otherwise the fiducial case. Also marked are the minimum 
mass for the presence of a shock-heated medium (solid, red) and 
the maximum mass for natural quenching by shocked accretion 
(dashed, magenta). 



ctd is rather weak. For a constant-density gas core, 
ctg = 0, any reasonable value for the total density slope 
in the range ad = — 2 would provide enough gravi- 
tational potential gain for balancing the cooling rate in 
Mv ^ 1O"'^^M0 haloes. Naturally, an isothermal density 
cusp (ttd = 2) provides a higher potential gain than a 
flat density core (ad = 0). We also find (not plotted) 
that the sensitivity to how deep in the potential well 
the energy deposit occurs is also rather weak, less than 
a factor of ^ 2 in Mv as long as the energy deposit oc- 
curs inside r < 0.5i?v The sensitivity to the gas profile 
is stronger, through the cooling rate. For an NFW po- 
tential well, the gas core has to be ag ^ 0.5 for H/Q > 1 
in ^ 10^'^ M© haloes. Once the gas profile is too cuspy, 
ag ~ 1, only haloes above lO^^-^M© can be quenched. 
If ag is as steep as 1.4, the cooling rate becomes so high 
that it cannot be balanced by gravitational heating even 
in cluster haloes of ~ 10^^ Mq. 

Figure [5] shows the dependence of Mmin on the mass 
fractions in hot gas versus cold clumps, /g and /c. The 
total gas fraction imposes a constraint on the way these 
parameters can vary, fg + fc — const. This could be at 
most the universal baryonic fraction, fg + fc — 0.15. A 
more realistic estimate would be fg + fc — 0.1, after 
subtracting the mass fraction of baryons in stars and 
those removed by feedback processes. With a total gas 
fraction of 0.1, one needs /c//g > 0.5 for quenching 
Mv > lO^^M© haloes. A ratio of /c//g > 0.1 would be 
enough for quenching Mv > lO^'^ ^M© haloes. We also 
find (not plotted) that the effect of increasing today's 
metallicity from Z = 0.3 to Z ~ Zq is only a ~ 50% 
increase in Mmin- 
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Figure 4. Estimated minimum halo mass for gravitational 
quenching Mmin as a function of the slopes of the gas and to- 
tal density profiles, ag and a^. The fiducial case is marked (red 
square) . 
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Figure 5. Estimated minimum halo mass for quenching M^nin as 
a function of the mass fractions in hot gas and cold clumps, /g and 
fc- The fiducial case is marked (red square). The relevant range 
is limited to below a diagonal line, e.g., /g + /c = 0.1 (dashed) or 
/g + /c = 0.15 (solid). 



The robust result is that for a broad and sensible 
range of values of the model parameters, once today's 
halo is of Mv ~ IO^^Mq or higher, there is in prin- 
ciple enough total power in the accretion available for 
keeping the halo gas hot. We conclude that the clumpy 
accretion scenario passes the first feasibility test: it car- 



ries enough gravitational power to overcome the overall 
cooling rate in haloes that host groups and big ellipti- 
cal galaxies, and more power than necessary in clusters 
of galaxies. The key is of course to have enough of this 
energy deposited in the inner halo, where most of the 
cooling tends to occur. This is addressed next, leading 
to constraints on the range of masses for the clumps 
that can serve as heating agents. 



3 HEATING BY GAS CLUMPS 

The actual accretion is aspherical and clumpy. The in- 
teraction of the infalling cold, dense gas clumps with 
the hot medium present in > IO^^Mq haloes, via ram- 
pressure drag and weak shocks, transfers energy into the 
gas and keeps it hot. The effect of dynamical friction, 
which is dominant for the most massive subhaloes, will 
be discussed in ^ 



3.1 Ram-Pressure Drag on a Clump 

In our analytic estimates, we consider a population of 
dense, cold gas clumps of mass each, falling radi- 
ally at the virial velocity through a halo of mass 
containing hot-gas mass Mg = fgM^. We assume that 
the clumps keep an internal temperature of ~ lO^K, 
where the atomic cooling function drops sharply and 
where the suppression of further cooling may be as- 
sisted by photoionization due to an external UV flux. 
The clumps are assumed to be in pressure equilibrium 
with a hot medium at the virial temperature, Tg — T^. 
Since the pressure in each of the gas components is 
P — pkT/M, the pressure confinement implies that at 
any position within the halo the densities of the cold 
and hot phases are related by 

PgTg = PcTc . (32) 
The implied clump radius is 

'^c = -^Y^s "^c, (33) 

describing the clump shrinkage as it moves into the 
higher density regions of the inner halo. 

The amplitude of the drag force acting by the hot 
medium on a cloud is 



/drag - TTCdPgK^ 



(34) 



where Cd — 1, indep endent of whether the m otion is 
sonic or supersonic (|Landau fc Lifshit"3 Il959l chapter 
XIlj3 Thus, the deceleration of a cloud due to drag 
in terms of the characteristic gravitational acceleration 
KV^v is 



^ The value of for a spherical clump varies from 0.4 to unity 
and back to 0.8 in the sub-sonic, trans-sonic and super-sonic 
regimes respectively. 
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(35) 



This implies that less massive clumps suffer a stronger 
deceleration, which naturally gets stronger with increas- 
ing ambient gas density. In turn, the rate of drag work 
done by the whole clump population of a given total 

mass is also oc rric . 

For a numerical estimate in the cosmological con- 
text, we use the top-hat virial relation, eq. ([9]), to obtain 



where 



1/3 ' 



T4 



me 



(36) 



(37) 



/Lt ~ 0.89. The deceleration can then be written as 
/< 



drag 



= 0.79 m, 



'1/3 



Rv 



(38) 



The work done per unit clump mass as it falls from 
radius fi to rf is then 



Din,ri) 



' /drag 



where / is a dimensionless integral 
For our fiducial gas profile this is 



i/3dr' , 



[3Ao(C)]-i/3ln 



1 



1 + a;f 



_ r 



(39) 



(40) 



(41) 



With C = 9, we have Ao{C) ~ 1.0 and Pg ~ 30 at r^. 
Then for infall from to rg the integral is / ~ 1.12 
and D ~ 0.79m^^/^K2 

Hydrodynamical simulations indicate that most of 
the work done by the ram-pressure drag is indeed de- 
posited i n the ambient gas rather than in the dens e 
clumps (jMurrav fc LinI |2004 iMcCarthv et all I2007D . 
This implies that the clumps can serve as efficient heat- 
ing agents, partly through the generation of turbulence 
(see ^J7]). The challenge is to have the clumps heat the 
gas in the inner halo, where its density peaks and the 
cooling rate is high. The clumps have to penetrate to 
the inner halo before they are stopped by the drag or 
destroyed by fragmentation due to hydrodynamical in- 
stabilities. This imposes a lower bound on the clump 
mass. On the other hand, the drag work has to be effec- 
tive such that it overcomes the radiative energy losses. 
This imposes an upper bound on the clump mass for a 
given halo mass. 



the core with a velocity higher than the virial velocity, 
accelerated by gravity despite the drag. Then it brings 
in more energy than is needed for heating the clump 
itself to the virial temperature when it melts into the 
hot medium. Such a penetration is possible only if the 
drag work done on the clump while it falls from i?v to 
Tg at a velocity ~ does not exceed the corresponding 
potential gain, 



D{Rv,r,) < A0(i?v,r-s). 



(42) 



If this inequality is not obeyed, the clump rapid in- 
fall velocity must be slowed down by the drag before it 
reaches the inner core. Substituting D from eq. ([55)1 , the 
constraint for penetration in a Hubble time becomes 



3 r3 



mc > 4.9 X 10'' M© p-" I 



2.2 X lO^'MfT 



(43) 



The second equality, here and below, refers to our fidu- 
cial case of Mv = lO^^M© at 2 = 00 This estimate 
assumes that the clump retains its initial mass as it 
moves through the halo. A higher lower bound may be 
obtained when clump fragmentation is taken into ac- 
count. 

The clump ploughing through the hot medium suf- 
fers Kelvin-Helmholtz and Rayleigh- Taylor instabilities, 
which eventually make it break into smaller fragments. 
Simulations show that this happens roughly once the 
clump has ploughed through a distance d correspond- 
ing to a gas column of ma ss comparable to the clump 
mass (jMurrav fc LinI I2004D . This is equivalent to the 
intuitive requirement that the impulse exerted by the 
drag, fdiiigd/Vv, is comparable to the momentum of 
the clump, rricVv, making it capable of breaking the 
clump up. We note that the corresponding timescale 
for breakup can actually be larger than the timescale 
associated with the fastest lii iear growing mode of the 
Kelvin-Helmholtz instability (jMurrav et al.lll993t) . The 
fragmentation is slowed down by the confining pressure, 
similarly to the slow-down that occurs when the clump 
is self-gravitating. 

Thus, a clump of initial mass rric on an inward radial 
orbit from radius rj would fragment at radius rf which 
obeys 



Pg{r) nr^^ dr . (44) 



Using the expression for the drag work, eq. (j39|) . this 
translates to fragmentation at an rf that obeys 

D{ru n) =0.5 . (45) 

A comparison of eq. (jl5|) with eq. (jl^ implies that 
clumps that could have made it to the halo core 
[by obeying eq. (j43p ] may fragment before they have 
reached the core. For our fiducial gas profile, eq. (ji5|) 
becomes 



ln(l + xf) = ln(l + Xi) - 0.91 Ao(C)^/3 p-^ mg^^ , 



(46) 



3.2 Penetration: Minimum Clump Mass 

A necessary condition for a clump to be effective in 
heating the gas at the halo core is that it penetrates into 



^ In order to reach the core in a Hubble time rather than in 
a crossing time R^/Vv, the average infall velocity should be ~ 
0.2 Vv, reducing the drag by a factor 0.2^ and thus reducing the 
minimum clump mass by a factor of 0.2^/3 ~ 1/3. 
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Figure 6. Expected clump penetration for different initial masses 
in a fiducial halo of Mv = 10^"^ Mq at 2 = 0. The blue curves 
show clump mass rric at radius r/R^. The clump is assumed to 
move from the virial radius inward on a radial orbit at a constant 
velocity Vv, and to break up to A'f = 2 fragments at a time 
(symbols). The solid red curve at the top marks rriBEi below 
which the gaseous clump is stable against its own self gravity 
( i|3.4l l. The dashed green curves refer to a local balance between 
drag heating and radiative cooling ( i|3.6|l : d = q (long dash) and 
5.5d = g (short dash). 

We now assume that at rt the clump fragments into 
A'f new clumps of mass rric/Nf each, and allow them to 
penetrate further till they fragment again by the same 
rule. If iVf ^ 1, the fragments do not penetrate much 
beyond the first rf. In order to penetrate from i?v to 
Vs before fragmentation occurs, the clump has to be as 
massive as 

TOc > 4.0 X lO^Mo ^3 ~ 4 X lO^Mo . (47) 

If, on the other hand, N{ is of order a few, the 
penetration following each fragmentation could be sub- 
stantial. For example. Fig. [5] shows the penetration of 
clumps of different initial masses at the halo virial ra- 
dius, under the assumptions of fragmentation into two 
pieces at a time, A'f ~ 2, at radii that are obtained 
by a repetitive use of eq. (P^ . in our fiducial halo. We 
see that in this case clumps of initial mass significantly 
smaller than rric ~ 10^ Mq can never make it to the 
inner halo. Clumps of rric ~ IO^Mq make it into the in- 
ner halo once, while clumps that are significantly larger 
penetrate through and orbit about the center more than 
once before they are completely destroyed. Thus, as long 
as the clump does not completely disintegrate in its first 
fragmentation event, our estimate of the lower limit im- 
posed by penetration with fragmentation is roughly 

mc ^ IO^Mq . (48) 

It is worth noting that the effectiveness of clumps 
that lie significantly above the minimum mass and thus 
safely make it into the core is not affected by the frag- 



mentation, as it is determined by the depth of the po- 
tential well. For small clumps, however, the fragmen- 
tation enhances the drag deceleration and can there- 
fore prevent a clump that would have otherwise reached 
the core from actually doing so. This tightens the lower 
bound that were obtained before considering the frag- 
mentation. 

The lower limit on TOc scales strongly with a, mostly 
through ij,^ and partly through P (via C cx a). At z = 2, 
for Mv = 10^^ Mq, the minimum rric drops by a factor 
of ~ 40 compared to its value at z = 0. The scaling 
with halo mass is such that the minimum TOc at 2 = 
drops by a factor of ~ 9 when the mass grows to 
Mv = IO^^Mq. The decrease of fi with increasing z or 
Mv originates from the corresponding increase in Tv. 
This increases the confining pressure, which makes the 
clump shrink (at a given mass), thus reducing the drag 
force (eq. l35|) . and allowing better penetration with later 
fragmentation. 

3.3 Effective Heating: Maximum Clump Mass 

Clumps that are massive enough to make it into the 
halo core should not be too massive; otherwise they fail 
to dissipate enough of their energy in a Hubble time. 
We can now address the heating-to-cooling balance with 
the heating computed from the actual drag work rather 
than the available potential well. Assuming that each 
clump is performing its act for a Hubble time as it orbits 
back and forth in the halo, and noting that the Hubble 
time is always about 5.5 times longer than the virial 
crossing time i?v/K, we can estimate the drag work per 
unit mass as 5.5 D{Rv, rg) and thus require for effective 
drag that 

5.5D(i?v,r,)Afe>£;hcat. (49) 
This translates to 

i?(i?v,r,) > A0(i?v,r,) (^5.5-^^ , (50) 

where H/ Q is the global heating-to-cooling ratio as eval- 
uated in eq. pO|) and plotted in Fig. [2l The implied up- 
per bound is a straightforward modification of eq. (|43|) , 

rric < SAxlO"^ Mq (^^^ ~ 3.0x IO^Mq , (51) 

where the last equality is for our fiducial case where 
H/Q = 2. 

The upper limit of eq. ([51]) is likely to be too tight. 
First, massive clumps are expected to speed up beyond 
Vv due to the gravitational pull to the halo center, caus- 
ing a significant enhancement in the drag force. If the 
velocity becomes twice as high along parts of the orbit, 
say, the instantaneous drag force is four times larger, 
so the upper limit on mc may increase by an order of 
magnitude. Second, the successive fragmentation can 
increase the drag work by clumps that start more mas- 
sive than this limit. According to the condition for frag- 
mentation, eq. (|45p . the number of crossings A'cross of a 
distance i?v before the first fragmentation can be esti- 
mated by iVcross D{Rv, Ts) ~ 0.5 V^, namely 
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Across — 0.63 mg/^ /i ^ / 



(52) 



This implies that even a clump as massive as lO^M© in 
a fiducial halo of IO^^Mq would fragment before cross- 
ing the distance i?v three times, and the consequent 
enhanced deceleration by drag would give it a chance 
to deposit enough energy by drag work in a Hubble 
time. The actual, more relaxed upper limit on rric will 
be evaluated more accurately using simulations below. 

Note that the condition for effective drag in eq. (|49)) 
does not guarantee that the apocenter of the clump's or- 
bit decays into the core in a Hubble time. The condition 
for this is rather 5.5 Z? > A0, so the corresponding up- 
per limit is lower by (H/Q)^. This implies that clumps 
near the upper limit can deposit enough energy dur- 
ing their passages through the core without being fully 
confined to the core. 

The scaling of the upper limit on nic in eq. (|5ip with 
Mv and z is dominated by the dependence of [H/QY on 
these parameters (Fig. [2]) which is partly balanced by 
the opposite trends of P discussed in ^'3.'2\ Changing 
the mass to Mv = IO^^Mq at z = 0, the maximum rric 
is increased by a factor of ~ 6 x 10"^. This change is 
driven by the dependence of the heating rate per unit 

mass on cx Mv^^. Moving to z = 2 at a fixed halo 
mass, the maximum rUc is increased by a factor of ~ 13. 

3.4 Bonnor-Ebert Stability: Maximum Clumps 

For the clumps to be effective in heating the inner halo 
gas they have to remain gaseous as they travel through 
the halo. This means that their inner pressure should 
support them against collapse due to their own self 
gravity. 

A self-gravitating pressure-confined isothermal 
sphere, with an equation of state P = pc^ (7 = 1), 
is stable on ce it is less niassive than t he Bonnor-Ebert 
(BE) mass (lEbertlll955t lBonnoil[T956L eqs. 3.5,3.6), 



1.18 ct 



(53) 



Q3/2 pl/2 ■ 

Writing the speed of sound within the clump as Cg = 
(k/m)Tc^ and the pressure confinement by the ambient 
ideal gas (7 = 5/3) as P — pg fcTg/m, we obtain 



msE ^ 1.6 X lO^Afo Tg-'/" T| , (54) 

where p-27 = Pg/10~^^g cm^'^. In the cosmological con- 
text, using pg from eq. (|2T|) and eq. H]), we obtain 

p_27 = 2.76 X 10-2 pg/.05 -4-3. (55) 

Using eq. © for Tg, and p from eq. ([37]), we obtain that 
the upper limit for BE stability is 

mc ^ mBE ^ 6.3xlO«M0 p-^/^ /gf p"" ^ lO^M© .(56) 

For our fiducial profiles, with C = 9, we have pg ~ 0.24 
at i?v, Pg — 30 at and pg ~ 243 at r = 0. Thus, 
TiBE ~ IO^Mq in the outer halo and msE ~ IO^Mq in 
the inner halo. The Bonnor-Ebert mass is shown as a 
function of radius within our fiducial model in Fig. [51 

We thus find that the maximum mass for Bonnor- 
Ebert stability is comparable to the maximum mass for 



drag efficiency as (under-) estimated in eq. (|5ip . This is 
an interesting coincidence, which allows a broad mass 
range for the clumps that are capable of serving as heat- 
ing agents. 



3.5 Tidal Disruption: Maximum Clump Mass 

The clumps moving through the inner halo experience 
tidal forces that may disrupt them before they manage 
to heat the ambient gas. We assume that disruption 
happens when the tidal force acting on the surface of 
the clump becomes comparable to the force associated 
with the confining pressure. 

Th e tidal force inside a host halo can be estimated 
by (e.g. iDekel et al.l[200l 



Jt^r 5— rc rric ■ 



(57) 



where r — 2 ~ d\n M/d\nr. For an NFW profile it is 
T ;S 2 in the outskirts of the halo, t = 1 at rg, and t 
as r — *■ 0. The pressure force is estimated by 



fp-nrlP, 



(58) 



where rc is given by eq. (j33[) and P = pkT/m is the 
confining pressure by the ambient gas. Assuming that 
the ambient gas is isothermal at the virial temperature, 
the condition for stability against tidal disruption, ft < 
fp, is 



rc 2'' ^^rV^ir) 



(59) 



With rc from eq. p3p . and the virial temperature from 
eq. ([HI), the condition for tidal stability becomes 



mc < 7.5 X IOHIqTI^^ f.osM^^^'A^/^ 



Rl T3/2y3(r) 



(60) 



For our fiducial profiles, with C — 9, the r-dependent 
term in the second line is ~ 0.03 at rg, where Pg ~ 30, 
y ~ Vv, and T ~ 1. The clumps have to be less massive 
than 2.3 x IO^Mq in order to get undisrupted into the 
core of a IO^^Mq halo at z = 0. This estimate is likely 
to remain roughly valid at smaller radii throughout the 
core, where the increase in pg and the decrease in V and 
in r roughly balance the decrease due to r^. Thus, our 
crude condition for stability against tidal disruption in 
the halo core is 



mc < 2.3 X 10«A/o T^^' /.05 M^{^ A^^^ 



(61) 



The upper limit imposed by tides in more massive 
haloes is more relaxed. 

We note that tidal disruption imposes an upper 
limit comparable to those imposed by Bonnor-Ebert in- 
stability and by drag inefficiency, all roughly demanding 
for AU ^ 10"M(T, at z = 



TTlr € 10*M 



© 



(62) 
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3.6 Local Heating versus Cooling Rates 

Given an accretion rate of appropriate clumps into the 
halo, we can crudely estimate the heating rate by drag 
in a single passage through a given point within the 
halo, in comparison with the local cooling rate. This 
will provide a lower limit to the actual heating rate, 
which can be higher due to repeating passages of the 
orbiting clumps in a given position. 

Consider a radial accretion of clumps in a rate Mc- 
The work done by drag per unit time and unit mass of 
ambient gas of density pg in a shell of radius r is 



47rr^ V Pg 



(63) 



The clump velocity v drops out, and using eq. psp we 
obtain 



J n T/2 -^^c Rl —2/3 -1/3 

d = 0.26cdV^ — H—p '-'niQ 



(64) 



Inserting the average cosmological accretion rate from 
eq. (fT7|) . the heating rate becomes 



/g 



1-1. (65) 



Dividing d from eq. by the cooling rate from 
eq. we obtain 



If 
q r 



Zc,3.9^Pg'/'m-'/' 



J8/3 f 

XYV_23^^^13 'Z225 ~r •' ■^'^ 4 • 
" Jg 



(66) 



The local heating rate balances the cooling rate, d/q 
1, for 



n6 

6.1 X IO^Mq -fpgir) 



8 J- 3 

-23 ""13 Y •'•05 ^4 • 



(67) 



If each clump passes 5.5 times in each radius over a Hub- 
ble time, the balance condition becomes 5.5 d/q = 1, 
and the corresponding nic becomes 5.5^ times higher. 
These values of mdr) are plotted in Fig. [5] for our fidu- 
cial case. We learn that the drag heating by clumps less 
massive than ^ IO^Mq can overcome the cooling in 
the inner halo in a single passage, while more massive 
clumps require multiple passages in the same location. 
The actual heating rate by such clumps will be evalu- 
ated by simulations below. 



3.7 Gas clumps in Dark-Matter Subhaloes 

3. 7. 1 Effect on the Clump Mass Estimates 

So far we have considered self-gravitating gas clumps of 
gas mass mc. As will be discussed below, these clumps 
may initially be embedded in dark-matter subhaloes. 



such that the total effective initial clump mass is > 
TTic. As long as the dark subhalo is attached to the gas 
clump, we can crudely evaluate its qualitative effects on 
our clump mass estimates. 

In terms of the overall energetics considered in ^ 
the total gravitational energy associated with the ac- 
cretion of the clumps into the center of the potential 
well could become larger by a factor that may be a 
significant fraction of m^/mc. The potential for over- 
coming the overall cooling is therefore higher, namely 
the threshold mass for quenching could become as small 
as ~ IQi^M©. 

The clump fragmentation as estimated in ^'6.2\ 
is driven by the hydrodynamical instabilities at the 
boundaries of the dense gas component with the am- 
bient hot gas, which may or may not be affected by the 
presence of dark-matter subhaloes. The subhalo poten- 
tial well is likely to slow down the fragmentation pro- 
cess, which may allow even smaller clumps to penetrate 
to the inner halo. The estimated minimum clump mass 
of ~ IO^Mq may therefore remain the same or become 
somewhat smaller. 

On the other hand, once the drag force remains the 
same while the clumps carry larger inertia, their orbital 
decay (! j3.3p should be less efficient, and the associated 
maximum mass for effective heating of the inner halo 
gas by drag would possibly decrease by a significant 
fraction of m^/mc- The potential for overcoming the 
cooling by drag at the first passage would be weak- 
ened. Thus, an effective heating by relatively massive 
gas clumps in subhaloes would rely on a quick separa- 
tion between the gas and dark-matter components. 

3.7.2 Separation of Gas and Dark-Matter Subhalo 

The ram-pressure drag is acting on the gas component 
only. Once the drag is stronger than the gravitational 
force tying the gas clump to its dark-matter subhalo, the 
two components detach from each other and the treat- 
ment of the gas clump as isolated is justified. Assuming 
an NFW density profile inside the subhalo, the maxi- 
mum restoring gravitational force per unit mass, /r, 
is obtained somewhat below its inner radius, r ^v/c, 
where is the subhalo virial radius and c is its effective 
concentration parameter. The condition for separation 
then becomes 



r-v/c 



< 



/dr 



(68) 



where Wv is on the order of the virial velocity of the 
subhalo. Expressing /drag in terms of /R^ of the host 
halo, eq. (|38p . and using the standard relation between 
virial quantities, {V^ / R^) / {v"^ / r^) — {M^/m^Y/^ , we 
obtain separation for 



,(r)VV^/^Mif 4^-10«Mo,(69) 



where cio = c/10 and /o.i = (mc/mv)/0.1. The second 
equality is for our fiducial case at Ts, where pg ~ 30, and 
assuming cio ~ /o.i — 1. This means that all clumps of 
gas mass rric ^ 10® are expected to get separated 
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from their subhaloes before or during their first pas- 
sage through the inner halo. In more massive haloes, 
the separation is even more effective. 



4 SIMULATIONS OF GAS CLUMPS 
4.1 Simulations: Method 

The analytic estimates of ^J3]are based on several crude 
approximations. In particular, the clumps were assumed 
to move inward radially and in a constant velocity. On 
the other hand, we learned that the orbital decay of 
the 1O^~^M0 clumps into the halo core is important 
for their ability to effectively heat the ambient gas. We 
therefore perform simulations of heating by clumps as 
they move on general orbits that properly respond to 
the gravitational potential well and to the gaseous drag 
and dynamical friction. The simulations include clump 
fragmentation and possible re-merging, and they allow 
imposing Bonnor-Ebert stability when desired. 



4.1.1 Halo 

The spherical halo profiles, including dark matter and 
gas in hydrostatic equilibrium, are adopted from ^2.11 
In our current simulations, these profiles are assumed 
to remain static as the cold clumps fall through the 
halo, ignoring for now the structural response of the 
ambient gas to the heating by clumps. Our fiducial halo 
has Mv = W^Mq and /g = 0.05. The total density 
profile is NFW (ad = 1) and the gas has a flat density 
core (ttg = 0). 



4.1.2 Initial Clumps 

Each clump of initial mass rUc starts at the virial radius 
of the host halo with an initial velocity of amplitude v 
and a direction set such that each of the two tangen- 
tial components equals ^fp times the radial component 
(i.e., /3 = and /3 ^ 1 for radial and circular orbits 
respectively). In our fiducial case for the simulations, 
we start with v — V^, and (3 — 0.23, corresponding to 
a typical subhalo orbit in dissipationless N-body sim- 
ulat ions, where the rat io of pericenter to apocenter is 
1:6 (|Ghigna et al.1119981 ). The clump temperature is as- 
sumed to be fixed at lO'^K until it disintegrates. 

In order to avoid artifacts that may be associated 
with a specific choice of the initial conditions, we sim- 
ulate 4000 clumps with slightly different initial condi- 
tions and average the results. The initial values of mc, v 
and (3 are drawn at random from normal distributions 
about the chosen mean values, with standard deviations 
Om = 0.5 mc, CTt, = 0.4 w and = 0.6 respectively, with 
the additional constraint that the parameters are all 
non-negative. The accretion rate is assumed to be given 
by eq. PT|) . with a fraction of mass in clumps in the 
fiducial case /c = 0.05, such that the total gas fraction 
is /g + /c = 0.1. 



4.I.S Forces 

Each clump is moved as a test particle by integrating its 
equation of motion in time through three-dimensional 
space. The forces that act on each clump are the grav- 
itational force exerted by the host-halo potential well, 
the drag force due to the ambient gas, eq. and the 
dynamical friction acting on the clump by the dark- 
matter and gas of the host halo (" OstrikeilllQQgf ). The 
clump radius entering the drag force is varied according 
to the pressure balance with the ambient gas, eq. ([55)1 . 
assuming that the temperature and density are uniform 
within the clump. 



4.1.4 Fragmentation and Mergers 

Fragmentation is assumed to occur once the clump has 
ploughed through ambient gas mass equal to its own 
mass. This is computed by performing an integral sim- 
ilar to eq. (|44| but along the actual orbit of the clump. 
At this point the clump breaks into iVf equal fragments, 
which are assumed to conserve energy and momentum 
and continue on the same orbit side by side. 

In some test runs, the clumps are assumed to merge 
when they collide with each other. These mergers could 
make some difference in the dense inner halo, where the 
collision rate could be high and the mergers may bring 
the fragmentation process to saturation. The mean time 
between collisions is estimated by t = (rtaw)"^, where 
n is the number density of clumps at that radius, a = 
irr'^ , and v is the clump velocity, all averages over many 
clumps on different orbits. Since the density of clumps is 
affected by the mergers themselves, it is determined by a 
numerically stable iterative procedure, where the clump 
number density profile is the average of the final profiles 
from the previous two iterations. The actual time for a 
collision event of a clump, dt, is drawn at random from 
an exponential probability distribution P{dt) oc e"''*/'^. 
Once a clump goes through a collision, the clump mass 
is increased by an amount equivalent to the average 
clump mass at that radius. The collision is assumed to 
be totally inelastic, where momentum is conserved and 
kinetic energy is lost. A certain fraction of it, /m, is 
assumed to be deposited in the ambient gas, and the 
rest is lost to radiation. 



4.1.5 End of Clumps 

We stop simulating the evolution of a clump once it ei- 
ther becomes bound within the innermost 0.002it!v of 
the host halo or it has been disintegrated to fragments 
smaller than 10~^ of its initial mass. The remaining ki- 
netic energy of the clump, minus the energy necessary to 
heat its gas to Tv, is added to the ambient gas at that 
radius. When the clump survives for a period longer 
than 10 Gyr, or when the clump exceeds the Bonnor- 
Ebert mass (if this feature is turned on), the clump is 
removed from the simulation without adding any fur- 
ther energy to the medium. 
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Figure 7. Time evolution of clump properties from a simulation 
of a single clump in our fiducial case. Shown are the clump radius 
in the halo r/Rv (solid, red), the clump velocity v/Vv (dashed, 
green) , the clump mass rric in units of its initial value at Rv (solid 
steps, blue), the fraction of the energy deposited in the hot gas 
(dot-dashed, cyan), and the clump mass relative to the critical 
mass for stability rriBE (dotted, magenta). 

4.1.6 Gas in Subhaloes 

When testing the case of gas clumps embedded in dark- 
matter subhaloes, we assume complete separation be- 
tween the two components once eq. is satisfied. 

4.1.7 Fiducial Case 

The fiducial case in the results presented below consists 
of gas clumps with masses about rric = IO^Mq in a 
halo of Mv = 10^^ Mq. The gas fractions in the hot 
medium and in the cold clumps are /g = /c = 0.05, 
with metallicity Z = 0.3. Clumps start at the virial 
radius with a mean velocity of w = Vv and /? — 0.23. 
The fragmentation is into N{ = 2 fragments, mergers 
are ignored, and BE stability is not imposed. 

4.2 Simulations: Individual Clumps 

We start by studying the evolution of a single clump 
along its orbit. Figures [71 to [TT] show the time evolution 
of the following clump properties: its radius within the 
halo, its velocity, its mass, the fractional energy lost 
by the clump (and deposited in the ambient gas), and 
the clump mass relative to the Bonnor-Ebert mass for 
stability. Figure [7] shows the fiducial case, and in the 
following figures we vary one parameter at a time about 
the fiducial case. Although the different cases may show 
different behaviors, we will see in the next section that 
in many cases they are as effective as cooling agents. 

In the fiducial case, the 10'' M© clump and its frag- 
ments move along an orbit through a sequence of peri- 
centers and apocenters. The orbit decays because of the 
ram-pressure drag, and it enters the 0.1i?v core by i ~ 4 
Gyr. The clump fragments for the first time at the vicin- 
ity of the first pericenter, near r ~ 0.2i?v at t ~ 2 Gyr, 
and it practically disintegrates while it is well inside the 
inner core by '--^ 6 Gyr. The drag is at maximum when 
the velocity and density are at maximum, near pericen- 
ter, and this is where the energy loss rate is maximal. 
By the time the orbit has decayed to r ~ 0.1i?v, the 
clump has deposited in the gas about one half of the 
total energy that it has potentially available, while the 
other half is deposited inside 0.1i?v by i ~ 6 Gyr. 
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Figure 8. Evolution of clump properties: dependence on frag- 
mentation scenario. No fragmentation (A^f = 1, top) versus strong 
fragmentation (A'^f = 10, bottom) to be compared to the fiducial 
case of A^f = 2 shown in Fig. [7| 

Note that the velocity remains for many orbits on 
the order of V^, to within a factor of 2, justifying the 
crude approximation made in ^ However, the factor 
of 2 increase in velocity during the first infall, cor- 
responding to a factor of 4 in drag force, makes the 
lO^M© clump more effective than estimated from Fig.[Sl 
Eq. ((5T|) is therefore an underestimate of the actual 
maximum clump mass for effective heating in a Hub- 
ble time, as suspected. 

We also see in Fig. [7] that the clump mass is signif- 
icantly below tobe at all radii, meaning that it remains 
gaseous and stable against collapse and star formation 
under its own gravity. 

Figure [5] explores the dependence on the fragmenta- 
tion recipe. In the top panel the fragmentation is turned 
off, A^f = 1, namely the clump mass remains constant 
as it moves along its orbit through a sequence of peri- 
centers and apocenters. The orbital decay into 0.1i?v is 
only slightly slower than for A^f = 2, but it does be- 
come less efficient well inside the core, reaching 0.01i?v 
at twice the time. Correspondingly, the total energies 
used for heating the core in the A^f = 1 and N[ — 2 
cases are comparable. The bottom panel presents a case 
with more efficient fragmentation, A^f — 10. The decay 
into the inner core is faster, but the clump completely 
disintegrates before it manages to make it into the very 
center. 

By simulating cases in which the clumps are as- 
sumed to merge when they collide, we find (not plotted) 
that the effect is negligible when the fraction of the en- 
ergy deposited in the hot medium is /m ~ 0.5. Even 
in the extreme case where the colliding clumps are as- 
sumed to disappear without depositing any energy into 
the hot medium, the effect of mergers is rather small. 

Figure [5| shows that the clump mass is significantly 
below Tobe at all radii for any fragmentation scenario, 
with the BE stability increasing with increasing A'f . We 
also find that the stability is slightly decreasing with 
more efficient merging. 

Figure [5| explores the effect of orbit eccentricity. As 
a function of decreasing /3, from circular (/3 2> 1) to 
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Figure 9. Evolution of clump properties: dependence on orbit 
eccentricity. An almost circular orbit (/3 = 100, top) versus a 
radial orbit (/3 = 0, bottom) to be compared to the fiducial case 
/3 = 0.23 shown in Fig. [3 
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Figure 11. Evolution of clump properties: dependence on halo 
mass. A galactic halo (A/v = IO^'^Mq, top) versus a cluster halo 
(Mv = 10"*^^ Mq, bottom), to be compared to the fiducial case of 
Mv = lO^'-' Mq shown in Fig. [7] 
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Figure 10. Evolution of clump properties: dependence on clump 
mass, from rric = 10* Mq (top) to rric = Mq (bottom), to be 
compared to the fiducial case of rric = IO^Mq shown in Fig. [7] 

radial {(3 = 0) orbits, the orbital decay rate and tlie 
fragmentation rate are both faster. This results in a 
higher heating rate in the core, but only by a factor less 
than three from circular to radial orbits. While for any 
eccentricity the clump mass is below mBE, it is getting 
close, mc/mBE ~ 0.5, in the case of a radial orbit near 
the first passage through the center. When starting in 
a radial orbit at rest, instead of with an infall velocity 



V = Vv, the first entry to the core happens later, but 
it has no significant effect on the evolution once inside 
the core (not plotted). 

Figure [To] shows the dependence on clump mass for 
Mv = IO^^Mq. At rric IO'^Mq, the clump disinte- 
grates well outside the core and deposits all its energy 
there. Such clumps cannot serve as heating agents for 
the core. Clumps in the range rric ^ IO^^'^Mq are all 
properly fragmenting, penetrating to the core and de- 
positing their energy there. The time it takes for the 
orbit to decay is longer with increasing clump mass, 
but all the clumps in this mass range deposit most of 
the available gravitational energy in time, and they are 
BE stable at all times. At rric ~ 10^ Mq, the orbit decay 
is slow so the clump barely makes it into the core in a 
Hubble time. Thus, only a fraction of the energy avail- 
able in the potential well is deposited in the core during 
the multiple passages of the clump there. The total en- 
ergy deposited, and its balance with the cooling rate, 
will be clarified in the next section when we consider 
a whole population of clumps given an accretion rate. 
For an initial clump of 10^A/q, the clump mass is only 
slightly below tobe when it enters the core. The orbit of 
an rric IO^Mq clump does not decay significantly in 
a Hubble time, so only about 10% of the energy avail- 
able in the potential well is deposited. This may lead to 
insufhcient total energy deposit for balancing the cool- 
ing. For an initial clump of IO^Mq, the clump mass is 
always above tobe- 

Figure [TT] shows the dependence on halo virial mass 
for rric = IO^Mq. The times for fragmentation and or- 
bital decay are increasing with halo mass. Still, the orbit 
does decay to inside the core in a Hubble time even for 
Mv ~ lO^^M©, and the vast majority of the available 
gravitational energy is deposited in time. 

Figure [12] shows the average properties of the clump 
distribution as a function of radius in steady state from 
a simulation of our fiducial case. The mass density in 
clumps is comparable to that of the ambient gas near 
r ~ Q.li?v. It is ~ 4 times lower at the virial radius, and 
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Figure 12. Average log properties of the clump distribution as a 
function of radius in steady state from a simulation of our fiducial 
case. Shown are the clump mass rric in units of its initial mass 
(dot-dash, magenta), the density in clumps pc compared to the 
ambient gas density pg (short dash, blue), and the gas density 
pg versus the mean gas density pg (long dash, green). The ratio 
of energy in turbulence and in thermal energy (horizontal curves, 
red) is shown for = lO^^, IQl^, IO^^Mq C ^gH . 




Figure 13. Heating rate (red) versus cooling rate (blue) at radius 
r in our simulated fiducial case. The accretion rate is the cosmo- 
logical average for 10^^ Mq haloes, of which a fraction fc = 0.05 
is assumed to be in ^ IO^Mq gas clumps. The heating beats the 
cooling everywhere, quite uniformly, by a factor of ~ 3 or more. 



it becomes an order of magnitude higher at the inner 
core. 



4.3 Simulations: Heating/Cooling Rate 

Our main results from the current simulations are ex- 
pressed in terms of the gravitational heating rate ver- 
sus the radiative cooling rate at every radius within the 
halo. The accretion rate of clumps is assumed to be a 
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Figure 14. Ratio of heating to cooling rates (H/Q) at different 
radii in our simulations. The fiducial case (solid red) is compared 
to other cases where one parameter is varied at a time. Here we 
test the dependence on the fragmentation scenario (top), the orbit 
(middle), and the presence of a dark-matter subhalo (bottom). 
The results are rather insensitive to these features of the model. 



fraction — 0.05 of the total average cosmological ac- 
cretion rate as approximated by eq. (|17p. with a fraction 
/g — 0.05 in hot gas. 

Figure [T2] shows the actual heating and cooling 
rates for the fiducial case, clumps about nic — 10'' 
in a halo of = 10^'^ Mq. We see that the heating 
overwhelms the cooling by a factor of three or more at 
all radii. 

Figure [14] tests the sensitivity of the heating-to- 
cooling (H/Q) ratio to the fragmentation recipe, the 
orbits, and the presence of a dark-matter subhalo. With 
no fragmentation, or fragmentation to only a few frag- 
ments at a time, the H/Q ratio is 3 or more everywhere. 
With more efficient fragmentation, N{ = 10, the heat- 
ing is almost the same outside 0.1i?v, it is still winning 
outside 0.06i?v, but is barely sufficient inside 0.06i?v 
(short dash, blue). When including the effect of merg- 
ers between clumps in a way that is most unfavorable 
to the cause of heating, namely when they are assumed 
to disappear with no additional energy deposited in the 
medium, H/Q is still of order unity even inside the core 
(dot-dash, magenta). The orbit dependence is weak — 
only a factor of two between radial and circular orbits. 
Even clumps that start on circular orbits end up beating 
the cooling by a factor of order two or more everywhere. 
The inclusion of dark-matter subhaloes makes only a 
negligible change for nic = l(f Mq. The effect is weak 
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Figure 15. H/Q ratio as a function of radius in our simula- 
tions as in Fig. 1141 here exploring the dependence on clump mass 
mc (top), halo mass (middle), and clump fraction /c with 
fs ^ f<^ ~ 0.1 (bottom). Each of these parameters is making 
a significant difference. The heating balances the cooling once 
mc ~ Mv > 6 X IQI^Mq, and /c//g ^ 0.5. 



because of the separation of the gas and dark-matter 
components by ram-pressure near the first pericenter. 
We conclude that the results from the adopted fiducial 
case could be interpreted as fairly representative, with 
no great need to worry about the actual recipes adopted 
for the fragmentation and mergers, the orbits, or being 
embedded in dark-matter subhaloes. 

Figure [T51 explores the H/Q ratio at the different 
radii as a function of clump mass, halo mass, or clump 
and hot gas mass fraction. The top panel confirms our 
estimate from that in a halo of IO^^Mq the effec- 
tive clumps are limited to the range IO^-IO^Mq. While 
clumps of ~ IO^Mq do it in a single passage, clumps 
in the range lO^-lO^M© provide the necessary heating 
via repeating passages at a given radius within the core. 
Nevertheless, the H/Q ratio for clumps in the effective 
range are very similar, rather insensitive to the actual 
clump mass. Note that TOc ^ IO^Mq clumps deposit 
enough energy in the core even though their orbits do 
not manage to decay into well inside the core in a Hub- 
ble time. Despite the fact that each of these clumps de- 
posits only a fraction of the energy available in the po- 
tential well, the drag they exert collectively while pass- 
ing repeatedly through the core transfers enough energy 
to balance the cooling rate, as argued in ^'d.'i\ Clumps 
of mc ^ IO^Mq fail to overcome the cooling in the inner 
halo; their weak deceleration by drag is not enough for a 




Figure 16. Ratio of gravitational heating versus cooling rate 
H/Q in the inner halo, r < O.lRv, as a function of halo mass for 
different clump masses, based on our simplified simulations. 



significant orbital decay, and the energy deposited dur- 
ing the transient passages through the core is below the 
cooling rate. Clumps of mc ^ 1Q'^^^MQ also fail; they 
cannot penetrate to the core because the drag they suf- 
fer in the outer halo and the associated fragmentation 
are too effective. 

The middle panel of Fig. [15] shows the depen- 
dence on halo mass for mc ~ lO^M© clumps. Recall 
that once the halo is less massive than the threshold 
for virial shock heating, ~ IO^^Mq, there is no two- 
phase medium, and the heating mechanism addressed 
here cannot work. The H/Q ratio is above unity at 
all radii for haloes of ~ 7 x IQ^'^Mq and above, 
as estimated in ^ As seen before, H/Q > 2 every- 
where for Mv ~ IO^^Mq, relevant to big ellipticals or 
small groups. It rises to H/Q ^ 100 for cluster masses, 
Mv ^ 1O^^M0. Note that the heating efficiency in the 
inner core becomes roughly the same for My = IQ^^Mq 
and IQ^^Mq haloes. This is because of the weaker de- 
celeration by drag and the slower fragmentation of the 
~ l{f Mq clumps in the former, where Tv is higher. 

The bottom panel of Fig. [15] explores the depen- 
dence on the fractions of clump mass and hot gas 
mass, /c and /g, assuming that the total gas fraction is 
/g + /c = 0.1. We see that the H/Q ratio is above unity 
everywhere as long as /c > 0.03, or crudely /c//g > 0.5. 

Figure [16] summarizes our results by displaying 
the total heating to cooling ratio in the inner halo, 
r < 0.1i?v, as a function of halo mass and for different 
clump masses. Here, unlike in previous figures, we also 
show cases that deviate from the fiducial case by the 
values of the two mass parameters. As seen before, grav- 
itational heating by clumps with initial masses in the 
range mc ~ 1O^~^M0 overcomes the cooling in the cores 
of haloes of Mv ^ 10^'^ Mq and above. We now learn, 
for example, that clumps of mc ~ IO^A/q could in prin- 
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Figure 17. Summary of constraints on the clump mass and halo 
mass for effective gravitational heating versus cooling rate H/Q in 
the inner halo, r < 0.1i?v, based on our simulations. The vertical 
line (transition from blue to yellow shading) marks the halo mass 
threshold for s hock heating above which a hot medium exists 
llBirnboim fc D ckcl 2003 ; JDckcl & Birnboim 2006). The dashed 
contours mark H/Q = 1,3,10,30, with the procedure for elim- 
inating clumps once they are Bonnor-Ebert unstable turned off 
(black) or on (red). The less constraining upper limit due to tidal 
disruption, eq. I I61II . is shown (dot-dot-dash, red). A very crude 
lower limit by heat conduction on clump formation inside haloes 
is indicated (dotted green). The permitted range where H/Q > 1 
is roughly Mv ^ 7 X lO^^M© and 10^ ^ mc ^ W^Mq. 
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Figure 18. Same as Fig. 1171 but at z = 2. The mass range 
for effective gravitational heating is broader, both for the halo 
mass and the clump mass: Mv ^ 2 X IO'^^Mq and 10* ^ rric ^ 
lO'^-^M©. 



ciple be effective in heating Mv > 10 A/q haloes. This 
is true, however, only when ignoring the upper limit im- 
posed on rric by Bonnor-Ebert instability. We see that 
these massive clumps do deposit some energy in the 
cores of haloes below IO^'^Mq, during their few quick 
passages through the core, but this is not enough for 
balancing the cooling rate there. At the small-mass end, 
it is interesting to see that clumps of W^'^Mq, which fail 
to heat the cores in haloes below do manage 

to heat the cores of rich cluster haloes ~ 10"'^^Mq. This 
is because the strong confining pressure due to the high 
virial temperature compresses the clumps into smaller 
sizes with reduced surface area. This weakens the de- 
celeration by drag and the corresponding fragmentation 
efficiency, and thus permits better penetration of these 
small clumps into the massive-halo cores. We also find 
that no clumps can overcome the cooling in the cores of 
Mv < 6 X 1O^^M0 haloes, despite the presence of shock- 
heated gas in halos above IO^^Mq. This is because the 
accretion rate into such haloes does not carry enough 
total power for overcoming the cooling rate, ^ 

Figure [T71 summarizes the constraints from the sim- 
ulations in the Mv-TOc plane, otherwise adopting the pa- 
rameters of the fiducial case. The contours mark equal 
H/Q values as integrated inside the 0.1i?v core. The 
lower limit on the halo mass for H/Q > 1 is robust 



near M,„i„ ~ 6 x 10 ^^M©, consistent with the esti- 
mates of Sj2] based on the global energy balance, and 
rather insensitive to the clump mass. This is because, 
in the range 10^ < nic < lO^M©, the clumps manage to 
share with the ambient gas most of the energy available 
to them when falling deep into the potential well. In 
Mv ~ IO-'^^Mq haloes, the requirement H/Q > 1 puts 
the lower limit for effective clumps at rric — lO^M©, de- 
creasing to mc ~ 3 X lO'^M© in haloes of Mv ~ lO^^M©. 
The upper limit for effective drag is rric ~ lO^M© at 
Mv ~ IO^^Mq, rising steeply with halo mass. 

The lower, red contours at the top refer to equal 
H/Q values in simulations where clumps are elimi- 
nated once larger than the Bonnor-Ebert mass, eq. ((53)) . 
The BE H/Q = 1 contour provides the practical con- 
straint, requiring an upper limit of rric — 6 x lO^M© at 
Mv ~ lO^^M©, which is rising slowly to rric ^ lO^M© 
at Mv ~ lO^^M©. This is slightly different from the 
crude estimate of eq. (j56p . shown as a line with a small 
negative slope, because the simulations involve a distri- 
bution of masses and the Bonnor-Ebert instability does 
not necessarily occur exactly at r^. As expected, the 
upper limit imposed by tidal disruption is somewhat 
weaker than the constraints imposed by H/Q > 1. 

Figure [15] summarizes the constraints in the Mv- 
rric plane at z = 2. The minimum halo mass for gravi- 
tational quenching is now lower, Mmin ^ 2 X lO^^^/©, 



© 2002 RAS, MNRAS 000.[Tlf24l 



18 A. Dekel, Y. Birnboim 



similar to the estimate in The minimum clump mass 
for penetration is now smaller than it was at z = by 
an order of magnitude, as expected from the scaling 
with /i'^ oc in eq. (j48p . The maximum clump mass for 
Bonnor-Ebert stability is only slightly smaller than its 
value at z = 0, consistent with the scaling in eq. ([56]) . 

where the z dependence of /5g via C (x a almost bal- 
ances the z dependence via /i"^. The maximum mass for 
effective drag is now higher by an order of magnitude, 
as expected from the scaling with in eq. ([STjl . and 
the strong dependence of on C cx a in eq. p6|) . The 
simulations thus confirm the estimates that the gravita- 
tional heating is more effective at higher redshifts, with 
the mass range broader both for the halo mass and the 
clump mass. 



5 HEATING BY DYNAMICAL FRICTION 

Dynamical friction (DF) is another important channel 
of energy transfer between the accreting clumps and the 
host halo. Here, the whole bound subhalo mass, domi- 
nated by its dark-matter component, contributes to the 
heating. However, this is roughly balanced by the fact 
that only a small part of the energy is deposited in the 
ambient gas, the rest being spent on "heating" the host- 
halo dark matter. Being proportional to the square of 
the total satellite mass, the dynamical friction is signif- 
icant only for subhaloes more massive than a few per- 
cents of their host halo mass, say ttidf ~ 0.05Mv. As- 
suming as an upper limit on the DF effect that the whole 
potential gain of infall from the virial radius is deposited 
in the inner halo, as in [J2l the maximum possible effect 
of heating by DF can be estimated as in eq. p^ . but 
with the gas-clump accretion rate Mc = fcM^ replaced 
by /Dp/gAfv, where /df is the fraction of the total ac- 
creting mass in the form of haloes more massive than 
mDF- (We tentatively assume here that the energy is 
divided between the ambient gas and the dark matter 
in proportion to their masses, but in the simulation be- 
low we do take into account the different response of 
the gas to DF.) The maximum effect of DF heating can 
thus be crudely estimated by /df times the maximum 
effect of ram-pressure or shock heating by gas clumps. 

The factor / pp can be estim ated using the expres- 
sions derived by iNeistein et al.l ((2006) for the cosmo- 
logical accretion rate based on the EPS formalism. The 
average total accretion rate onto a halo of mass at 
time t is well approximated by 



1 dMv 
Mv dw 



< moF) ^ \/ - [^'(mDF) -a'(Afv)]"'/' .(71) 



1 dAfv 
Mv dw 



-[n\M,/q)-n\M^)\ 

IT 



1/2 



(70) 



with q ~ 2.2. Here a{M) is the rms linear density fluc- 
tuation encompassing mass M normalized by as, and 
w = 1.68/ D{t) with D(t) the growth rate of linear den- 
sity fluctuations in the assumed cosmology. In the rel- 
evant halo mass range, 1O^^~^^M0, this approximation 
is slightly higher than the explicit estimate in eq. pT|). 
but only by less than 15%. The part of this accretion 
that comes in haloes of mass below some mass mpF is 
approximated by 



For a standard ACDM cosmology, halo masses in 
the range relevant for galaxies and clusters, and mpp 
between 1% and 10% of Afv, we obtain quite robustly 
/df ~ 0.5. This implies that the heating by DF can con- 
tribute a significant fraction of the total gravitational 
heating. The maximum effect of DF in terms of heating 
versus cooling can therefore be described in Fig. [5] by a 
curve similar to the curve shown for the total heating 
but a factor of two lower, implying that a most effective 
DF heating may by itself balance the overall cooling at 
z = in haloes of JVU > 10^^ Mq. 

Based on such simple considerations and sim- 
plified simulations, the DF has indeed been ar- 
gued to have an important overall effect on slow- 
ing down the cooling in clusters, though perhaps not 
enough for prop erly balancing the cooling rate in 
the cluster cores ^ El-Zant et al.ll200l iKim et al.|[2005l : 
iKhochfar fc Ostriker ,2007i) . A similar conclusion can 
be obtained from the apparent failure of certain hy- 
drodynamical simulations of massiv e haloes to pre- 
vent cooling flows in all cases (e.e. iMotl erail l2004 
iBurns et alj|2007l : iFaltenbacher et al.ll2007f ). This is as 
long as they do not incorporate another heating source 
such as AGN feedback that i s practically put by hand 
(e.g. ISiiacki fc Springel l2006l and references therein). 
On the other hand , most of the cluster simulations by 
Nag ai et all (|2007D seem not to show substantial cooling 
flows (private communication with A. Kravtsov), which 
are somehow prevented by gravitational heating alone, 
without any additional energy source. Despite the fact 
that all these hydrodynamical simulations are expected 
to resolve the massive subhaloes that provide most of 
the DF work, the apparently conflicting results indicate 
that the dust has not settled yet on the actual role of 
DF heating. As far as heating by gas clumps, the cur- 
rent simulations clearly fail to resolve small enough gas 
clumps for the ram-pressure heating to play a major 
role. 

In order to complete our crude estimate for the pos- 
sible role of heating by DF, we use simplified simulations 
similar to those described in where the response of 
the ambient gas to dynamical- f riction is properly com- 
puted following lOstrike j (|l999f ). Since our simulations 
do not incorporate tidal mass loss from the clumps or 
the subhaloes that may be attached to them, they tend 
to overestimate the effect of DF, and should therefore 
be interpreted only as upper limits on the role of DF, 
much like the energetics estimates discussed above. We 
perform simulations similar to those described in 21 
but with the accreting subhalo masses drawn at ran- 
dom from the mass function derived from the accretion 
rates in eq. ([70]) and eq. ([7T|) . Figure [12] shows the re- 
sulting H/Q ratio in the 10% halo core as a function 
of halo mass, compared to the H/Q shown earlier in 
Fig. [Tni for our fiducial gas-clump model. We see that 
the DF heating can be as large as one third of the ram- 
pressure heating. Nevertheless, recall that the subhaloes 
providing the DF heating are automatically available 
in a standard cosmological accretion, while the origin 
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Figure 19. Maximum lieating by dynamical friction. Shown is 
the ratio of gravitational heating versus cooling rate H/Q in the 
inner halo, r < 0.1 -Rv, as a function of halo mass. The simulation 
method is as described in i|4.1l Our fiducial model for heating by 
cold gas clumps via ram pressure (red, solid curve) is compared 
to the heating by a cosmological distribution of dark-matter sub- 
haloes via dynamical friction (blue, dot-dashed curve). The esti- 
mates for the DF heating are upper limits because the weakening 
of the DF due to the tidal stripping of the subhaloes is ignored. 

of the population of gas clumps required for efficient 
ram-pressure heating is more speculative (i i6.2p . 

Figure [^D] compares the maximum DF heating by 
a single subhalo of mass mv (with the tidal striping ig- 
nored) to the ram-pressure heating by a gas clump of 
mass mc = /cWv. We see that DF heating is negligible 
for the clumps that are relevant for ram-pressure heat- 
ing, TOc = 10^ — IO^Mq, but it becomes the dominant 
effect for massive subhaloes that are on the order of 1% 
of their host halo mass or larger. 



6 DISCUSSION 

6.1 On the Physics of Dissipation 

The actual process of energy deposition by ram pressure 
should be investigated in m ore detail. This study can 
build upon the key resu lts of iMurrav fc LinI l|2004l ) and 
iMcCarthv et al.l ()2007D that most of the ram-pressure 
drag work is deposited in the ambient gas rather than 
in the dense clumps. When a cold clump moves sub- 
sonically through the hot gas, the energy is deposited 
as kinetic energy. The population of clumps moving on 
different orbits create turbulence on a scale comparable 
to the mean distance between clumps, which cascades 
down to the scale where viscous heating becomes effi- 
cient. If a non-negligible fraction of the energy is tem- 
porarily stored in the turbulence reservoir it can have 
interesting consequences on the efficiency of clumps as 
quenching agents. The turbulence increases the effective 
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Figure 20. A comparison of the heating by a single clump due 
to dynamical friction and ram-pressure drag, based on our simu- 
lations. The gas-clump mass involved in the ram pressure is mc. 
The total subhalo mass involved in the dynamical friction is as- 
sumed to be rriv = 20 rric . The ram pressure is dominant for small 
gas clumps, while the DF becomes more important for massive 
subhaloes. 

pressure on scales larger than the mean separation be- 
tween clumps, helping the halo gas expand in response 
to the clumpy accretion and hence slowing down its 
cooling rate. The turbulent component does not cool 
like the thermal component, and, since the clumps are 
smaller than the smallest eddies, the turbulence does 
not add to the pressure confining the clumps. 

Following iKolmogorovl (|1941h . in a steady state, 
with energy input rate per unit mass e into eddies 
of scale L, the energy stored in the turbulence per 
unit mass can be written as e turb ~ co(eL)^/'^, with 
cn - 2.1 (Landau fc Lifehitg ([l95l §32) and IPopoI 
(|2000l §6.5.2)]. The corresponding pressure is -r = 
(7— l)/9geturb- Using our simulations, we obtain at every 
radius the energy deposit rate e and the mean clump 
separation L, and compute the ratio of eturb to the ther- 
mal energy. This is plotted in Fig. [T^l We find that 
this ratio, at r = 0.1i?v say, ranges from ~ 20% to 
~ 5% when the halo mass ranges from 10^^ to lO^^M©, 
respectively. This crude estimate is on the same or- 
der of magn itude as result s of hydrodynamical simu- 
lations (Wis e fc Abeil l2007t [Kltenbacher et al. 2007j), 
perhaps indicating that clumpy accretion may play a 
non-negligible role in driving turbulence in these simu- 
lations. The non-negligible energy stored in turbulence 
at Mv ~ 10 ^^M© may reduce the lower limit for ef- 
fective quenching to below the Mmin ~ 7 x lO^^M© 
obtained without the turbulence reservoir. 

The fragmentation of the gravitating clumps due 
to hydrodynamical instabilities should also be studied 
in more detail , foll owing up on the pioneering work of 
iMurrav fc Lh] (|20b4. ). Simulating a single clump falling 
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under gravity through ambient hot gas may be feasi- 
ble immediately. Another challenge would be to inves- 
tigate how the system of clumps may develop a steady- 
state mass function above a minimum mass due to the 
competi ng effects of fragmen tation and coalescence (see 
hints in lLin fc Murravll2000l ). 

Our estimates of the cooling rate were limited so far 
to the radiative losses from the smooth hot medium, but 
in a following study one should also evaluate the pos- 
sible enhancements to the cooling rate associated with 
the clumps themselves and with the two-phase medium. 
One such cooling source is the gas stripp ed from the 
surface of the clump by ram pressure (e.g. iKlein et al.l 
I1994 for stripping when a clump is hit by a shock). 
This gas is expected to be at lower entropy than the 
ambient gas and therefore to radiate more efficiently. 
On the other hand, this energy loss may be limited by 
the coalescence of the s tripped gas and sma ll clumps 
back into bigger clumps (jLin fc Murravl[2000[ ) . Another 
potential source of energy loss from the hot medium is 
thermal conduction into the cold clumps, which strongly 
depends on the uncertain strength of th e conductivity 
(see estimates in lMaller fc Bullocdl2004l ). 

The Bonnor-Ebert analysis for the stability of 
the gas clumps against their own self gravity should 
be expanded to include thermal stability against the 
formation of a two-phase medium inside the clumps 
themselves, which is possible for particular values of 
metallicity (jSternberg et al.ll2002l . for Z - 0.1). This 
is potentially another route to star formation, which 
may weaken the effectiveness of the clumps as heating 
agents. 

6.2 On the Origin of Gas Clumps 

The origin of gas clumps in the desired mass range, 
containing a large enough fraction of the accreting gas, 
remains a key open issue for the feasibility of the model. 
At this point we only attempt a preliminary discussion 
of certain options. 

One possibility is that the gas clumps form em- 
bedded in the cosmological population of small dark- 
matter haloes that become subhaloes through a se- 
quence of minor mergers building up the larger halo. 
For the gas clumps to be in the proper mass range 
for effective ram-pressure, 10^~^Mq, the halo masses 
have to be in the range 10^~^Mq. As long as an ex- 
ternal ionizing flux is effective in keeping the gas at 
~ 10'*K, even in the vicinity of big haloes where the 
density is ~ 100 above the universal mean, gas cannot 
accumulate in haloes less massive than ~ 3 x lO^M© 
(|Gnc din 20001) and cannot survive evaporation from 
smaller haloes pBarkana fc Loebir200ll : lLoeb fc Barkanal 
[200lt iShaviv fc Deke]|l200ir Thus, proper clumps can 
exist inside subhaloes only if the ionization is ineffec- 
tive. This may be the case before the universe becomes 
reionized at z ^ 9, by which most ~ IO^Mq haloes are 
already in place (jMo fc Whitdl2003 . based on the Press- 
Schechter formalism). The p hotoionization may b ecome 
ineffective again after z ~ 2 (jBabul fc Reeslll992D . or m 
the likely event that the dense clumps become shielded 
from the ionizing flux. If the clumps are not ionized and 



the environment is still partly expanding, the gas can 
cool adiabatically, which may allow it to remain bound 
even in small haloes. 

We estimated fi i3.7.2p that even if the clumps come 
in as part of subhaloes, the dark component should not 
have a substantial effect on the role of the clumps as 
heating agents, since clumps in the relevant mass range 
would separate from their subhaloes by ram pressure 
before they enter the inner halo. 

Another, perhaps more relevant possibility is that 
the cold gas clumps fragment from the hot gas by ther- 
mal instability and establish a two-phase medium in 
pressure equilibrium. With the peak cooling rate at 
^ lO^K, this requires an ambient gas of > lO^K. Most 
of the gas is shock heated to such temperatures in haloes 
more massive than ~ 10^^ Mq, as well as in nearby 
collapsed pancakes an d filaments that feed such haloes 
l|Cen fc Ostrikeij|2006f ). Note that the clumps can serve 
as effective heating agents that balance the cooling in 
the halo core even if they form inside the halo. Once 
they fall into the halo center, the gain in gravitational 
potential is not much smaller than the energy provided 
by external clumps. We verified that falling in from rest 
does not change the results by much. 

Clumps can form by thermal instability as long as 
the cooling function A(T) is a decreasing function of T. 
This is valid for haloes with < lO^K (for ~solar 
metallicity), namely < IO^^Mq at z = 0. This 
mechanism cannot provide the clumps necessary for 
heating rich clusters of galaxies, unless the clumps form 
in the cooler filam ents outside the ha l o and then fall into 
the virial radius. iMaller fc Bullockl (|2004 MB04) pro- 
vide a detailed discussion of the formation and survival 
of such clumps, in the context of the formation of big 
disk galaxies in dark-matter haloes. They specifically 
address a lower limit to the clump mass due to ther- 
mal conduction, whose strength is parameterize d by /s, 
the fraction of the classical Spitzer conductivity (^pitzeg 
[1962). They estimate that conductivity would suppress 
the formation of small clumps and impose a lower limit 
of 

me ^ 4 X IOHIq mI^/'^ f,l(^ (72) 

(MB04, eq. 31). This impl ies that if the conductivity i s 
as high as fs ~ 0.2 (e.g. IZakamska fc Naravanll2003l ). 
it would limit the mass range of effective clumps that 
could form by thermal instability inside virialized haloes 
of IO^^Mq to mc 1O'^"*M0. On the other hand, if 
fs is significantly smaller, as expected in the presenc e 
of a uniform mag netic field (iChandran fc Cowlevlfl99l . 
the minimum mass imposed by conductivity would have 
a weaker effect on the allowed mass range for effective 
clumps. The lower limit for clumps that could form in- 
side the virial halo, with fs = 0.01 in eq. (|72p . is marked 
in Fig. 1171 The minimum mass imposed by conductive 
evaporation is smaller by an order of magnitude (MB04, 
eq. 36). Note that the limit of eq. (|72p is not valid for 
the formation of clumps in the pancakes and filaments 
at the vicinity of haloes. 

One might expect that limits imposed by forma- 
tion and survival arguments on the masses of cold gas 
clumps can lead to an enhancement of the gas fraction 
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in clumps that lie in the mass range that allows effective 
heating by drag. For instance, a fraction of the cold gas 
that has been prevented by conductivity or coalescence 
from being in clumps much smaller than ~ IO^Mq may 
find itself instead in more massive clumps, which can 
more effectively penetrate into the inner halo and heat 
it f ji3.2p . The gas in clumps above the stability thresh- 
old of TUBE ~ 10^ Mq, where stars can form, is likely to 
be re moved by supernova feedback after the first star- 
burst (iDekel fc Silklll986l : lDekel fc Wooll2003D . and thus 
be added to the reservoir of gas available for condensing 
into smaller, stable gas clumps that are more effective 
in drag-heating the gas. 



6.3 Where are the Clump Remnants? 

If nic ~ 5 X IO^Mq cold gas clumps are responsible 
for quenching in haloes of ^ 10^^ Mq, one expects 
a total mass in clumps of ~ 5 x IO^^Mq (/c/0.05) A/13 
acting over a Hubble time in each halo, and a small frac- 
tion of this mass in today's remnants of these clumps 
that are potentially detectable. This fraction should be 
significantly smaller than 20%, given that (a) the cross- 
ing time from the virial radius to the halo center is 
^ 0.2iHubbie, (b) the accretion rate is decreasing with 
time, and (c) clumps are destroyed on their way in. In 
order to quench a ~ 10^"^ Mq halo one needs about a 
third of the fiducial fraction of mass in clumps, namely 
fc ^ 0.2 (e.g. from Fig.[5]or Fig. [13]). The corresponding 
average accretion rate of clumps today at the virial ra- 
dius is Mc ~ lOAf0 yr^^, and it gets smaller at smaller 
radii within the halo as many clumps are destroyed on 
their way in. The actual accretion rate of clumps into 
the inner halo depends on the mass function of clumps. 
The flow of clumps toward the center may not be eas- 
ily detected as a cold flow because the clumps that heat 
the center are eventually destroyed, heated and blended 
with the hot medium. The flow rate of clumps toward 
the center can be signiflcantly lower than the cooling 
flow implied by the gas density and cooling time de- 
rived from X-ray observations as long as these clumps 
bring into the center a significant fraction of the energy 
gained by infall into the potential well. Indeed, in our 
simulations these clumps enter the central region with 
a velocity that is typically twice as high as the virial 
vel ocitv fFigs.lTlfTTI). 

iMaller fc Bullockl (|2004 §7) have argued that the 
observed population of High V elocity Clo uds (HVCs) 
in the halo of the Milky Wav (IOortill966D or the Lo- 
cal Group (|Blitj l2002t iMalonev fc PutmanI l2003l and 
references therein) is consistent in many ways with the 
required population of clumps. This includes the gas 
temperature of ~ lO^K indicated by a median FWHM 
line width of Aw ~ 25kms~^, the estimated typical 
cloud size and mass of ~ 5 x 1O^M0, the kinematics of 



the HVC population within the halo, and the estimated 
total number of a few thousands HVCs, implying a total 
mass of ~ IO^^Mq in a ^ IO^^Mq halo, in the ballpark 
of our model predictions. 

MB04 (§8) have also argued that the desired clumps 
may be detected at high redshift as the high column 



density absorption systems in quasar spectra, such as 
Lyman l imit and CIV syste ms (based on the model of 
IMo fc Miralda-Escudd[l996l ) . which are observed to re- 
side within the ex tended virial haloes of bright galaxies 
lichen et alllMl . 

The remnants of the desired clumps may be de- 
tected as a minor component of cold gas in other- 
wise quenched ellipticals, groups and clusters. A large 
fraction of today's ellipticals contain non-negligible 
amounts of cold gas, detected by radio observations 
as HI or CO, a t the level of IO^-IO^Mq per galaxy 
(e.g. ISage et all |2007[ ). A significant amount of HI 
seems to be present in the inter-galactic medium inside 
groups of galaxies (L. Blitz, private com munication) . 
Ionized, warm ga s detecte d in ellipticals (jLauer et all 
l2005HSarzi et al.i r2006: G raves et al.ll2007D . sometimes 
as LINERs, may be "the tip of the iceberg" indicat- 
ing a significant component of neutral gas (S. Faber, 
private communication). Finally, the clump remnants 
may be associated with the "cold cores" detected in 
about half of the X-ray clusters, in which the tempera- 
ture is slightly co oler than the virial temperature (e.g. 
IChen et al.ll2007l . and references therein). 



7 CONCLUSION 

We performed a preliminary feasibility study of a sim- 
ple mechanism — gravitational heating by accretion — 
for long-term quenching in dark-matter haloes above 
a threshold mass ~ 10^^ Mq. Such quenching is de- 
sired for explaining the existence of red fc dead el- 
liptical galaxies and the lack of cooling flows in clus- 
ters of galaxies. It is now common wisdom that the 
halo gas is flrst heated to the virial temperature by 
a global s hock once the halo grows above this thresh- 
old mass (iBirnboim fc Dekell l200i iKeres et al.l 120051 : 
IDekel fc Birnboimll2006r and references therein). While 
the interaction of the accreting gas with the rapidly ex- 
panding shock is e nough for long-term q uenching in ha- 
los of 1012-13Mq ijBirnboim et al.ll2007D . we addressed 
here the possibility that the long-term quenching main- 
tenance in more massive haloes is due to the gravita- 
tional energy of cosmological accretion, being delivered 
to the inner-halo hot gas by cold flows, partly clumpy, 
via ram-pressure drag and local shocks as well as dy- 
namical friction. A robust estimate reveals that with 
the average cosmological accretion rate onto haloes of 
^ 1O-^^M0, the gravitational energy power available by 
gas infall into the bottom of the potential well can bal- 
ance the overall radiative cooling losses. This is a nec- 
essary condition for effective heating, independent of 
the actual mechanism by which the energy is deposited 
in the inner hot gas. This estimate makes gravitational 
heating a viable competitor to the fashionable alterna- 
tive of AGN feedback, and motivates a further study. 

We then addressed specifically the case of clumpy 
accretion, in which cold gaseous clumps are pressure 
confined by a hot medium. We found that this is a 
feasible quenching mechanism in haloes of a few times 
lO^^AfQ and above provided that the gas clumps are 
in the mass range 10^~^Mq. Smaller clumps are slowed 
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down by the ram-pressure drag and they disintegrate 
due to hydrodynamical instabihties before reaching the 
inner halo. They do not reach the core with a sufficient 
excess of kinetic energy. In addition, their formation 
could be suppressed by heat conductivity. On the other 
hand, gas clumps that are too massive do not trans- 
fer enough energy by ram-pressure drag in a Hubble 
time. Independently, they may stop being gaseous be- 
cause they are unstable to collapse and star formation 
under their own self-gravity, and they are more sus- 
ceptible to tidal disruption in the inner halo. By sim- 
ulating the process with gas clumps of a proper mass, 
we confirmed that the heating rate can indeed balance 
the cooling rate throughout the whole halo, as long 
as Mv > Mniin ~ 7 X lO^^M©, the gas inner density 
cusp is not steeper than pg oc r~°-^, and the mass frac- 
tions in cold clumps and in the hot ambient gas are of 
the same order of magnitude. The effect is stronger at 
higher redshifts, thus possibly affecting the core struc- 
ture in a way that makes the maintenance easier also 
at later times. The cold flux into the center is expected 
not to introduce a cooling-flow problem by itself be- 
cause the accelerated clumps carry into the center more 
energy per unit mass than the thermal energy of the hot 
medium and because they dissolve into the hot medium. 
We conclude that the gravitational quenching scenario 
has passed successfully several non-trivial preliminary 
feasibility tests. This refers in particular to the version 
where the energy is dissipated via rather low-mass cold 
gas clumps. The more massive dark-matter satellites 
can have a substantial additional contribution via dy- 
namical friction. 

Many of the elements of the proposed scenario 
should be investigated in more detail via hydrodynam- 
ical simulations, where the challenge of properly re- 
solving ~ 1O^~^M0 gas clumps in ~ IQ^^Mq haloes 
is not trivial. While there have been preliminary at- 
tempts to re solve small clumps using SPH simulations 
(|Kaufmann e t al. 2006) , a proper treatment of the hy- 
drodynamical instabilities associated with the clumps 
will probably require the use of an adaptive Eulerian- 
grid simulation technique (jAgertz et al.ll2007l ). 

An immediate next step could be performed using 
hydrodyna mical simulations , start ing with the spheri- 
cal code of iBirnboim et al] (|2007D . While the current 
study has been restricted to clumps falling into a static 
halo, it would be desirable to incorporate the dynamical 
response of the halo gas to the over-heating in excess of 
the cooling rate. We expect the subsequent expansion 
of the gas to slow down the cooling and make the heat- 
ing even more effective. In order to match the observed 
X-ray gas in massive haloes, the system will have to 
relax into a steady-state configuration with a constant- 
entropy core. It is possible that turbulence and convec- 
tion may be the key stabilizing processes. 

In parallel, cosmological simulations that zoom in 
on massive haloes will allow an analysis of heating by 
the smoother, massive cold flows. These studies should 
reveal whether the observed entropy floor in cluster 
cores can indeed be reproduced by gravitational heat- 
ing. If so, it may alleviate the need for special scenar- 
ios based on either smooth accretion after pre-coUapse 



heating (jVoit k Ponmanll2003t [Borgani et aT1l2005D or 

AGN feedback as mechanisms for enhanced entropy 
production. 

Based on our preliminary feasibility tests we con- 
clude that gravitational heating, e.g., by clumpy accre- 
tion, is a viable scenario, complementary to AGN feed- 
back as a long-term quenching mechanism. They are 
both triggered by the shock heating of the halo gas 
once above the threshold mass of ~ lO^^M©. While 
the energy transfer from the small black-hole scales to 
the extended halo gas requires a non-trivial and yet un- 
known physical mechanism, the appeal of the gravita- 
tional heating scenario is in its simplicity and availabil- 
ity. In particular, this mechanism provides a natural ex- 
planation for the characteristic halo mass above which 
quenching is effective, and for the uniform heating over 
the whole halo. 

The indicated dual role of cold flows in galaxy for- 
mation is fascinating. In haloes below ~ lO^^M©, where 
the cooling time is shorter than the relevant dynam- 
ical time, they seem to serve as the main source for 
disc buildup and star formation. In the massive haloes, 
where the cooling rate is slow, a stable virial shock can 
be sustained and the potential well (cx V^) is sufficiently 
deep, the cold flows, clumpy or smooth, can become 
the destructive quenching agents. They keep the halo 
gas hot and thus allow the suppression of gas supply to 
the central galaxy and the removal of gas from satel- 
lite galaxies. The potentially important role played by 
cold flows, both in the buildup of galaxies and in their 
subsequent quenching, deserves further detailed inves- 
tigations. 
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